Chemical Equilibrium in Heavy Ion Collisions: 
Rapidity Dependence. 



F. Becattini 1 , J. Cleymans 2 

1 1, Universita di Firenze and INFN Sezione di Firenze 
Largo E. Fermi 2, 1-50125, Florence, Italy. 

2 ,UCT-CERN Research Centre and Department of Physics, University of Cape 
Town, ZA-7701 Rondebosch, South Africa 

E-mail: becattini@fi.infn.it Jean. CleymansOuct . ac . za 

Abstract. Particle yields in heavy ion collisions show an overwhelming evidence 
for chemical or relative chemical equilibrium at all beam energies. The rapidity 
dependence of the thermal parameters T and \ib can now be determined over a wide 
range of rapidities and show a systematic behavior towards an increase in [i^ away 
from mid-rapidity. 



1. Introduction 

Over the past decade, the analysis of particle multiplicities in heavy ion collions has 
shown overwhelming evidence for chemical equilibrium in the final state except for 
particles carrying strangeness which are mildly suppressed; however, their relative yields 
fulfill statistical equilibrium. A summary as of 2006, combining the results from many 
different groups [lj, is shown in Fig. [TJ Except for particle multiplicities at RHIC 
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Figure 1. Temperature vs. [Ib us determined from heavy ion collisions at different 
beam energies. The lower AGS points are based on a preliminary analysis of 4-7T data. 
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energies, all data in Fig. [T] use integrated particle yields, the very systematic change of 
thermal parameters over the full range of beam energies is one of the most impressive 
features of relativistic ion collisions to date. It is now possible to use the thermal 
model to make solid predictions for particle multiplicities at LHC energies [2j and to 
determine which beam energy will lead to the highest baryon density at freeze-out [3]. 
With chemical equilibrium for integrated yields thus firmly established, we focus on other 
properties, in particular, since the rapidity distributions of identified particles is now 
becoming available also at RHIC energies [7] , it is now possible to determine the rapidity 
dependence of thermal parameters. A first analysis was done by Stiles and Murray [I] 
for the data obtained by the BRAHMS collaboration at 200 GeV [61 [7J . This shows a 
clear dependence of the baryon chemical potential on rapidity in particular due to the 
changing p/p ratio. Subsequently it was shown by Rohrich [5] that the particle ratios at 
large rapidities are consistent with those measured at the SPS energies. This opens the 
possibility to compare measurements for e.g. the K + /n + ratio at high rapidities and 
check them with the corresponding values measured in the energy scan at the SPS, thus 
complementing the rapid variation of this ratio as a function of beam energy. A further 
analysis of the rapidity dependence was recently done in Ref. [8]. A detailed report of 
our results will be presented elsewhere [9]. 

2. Rapidity Distribution 

The general procedure is as follows: the rapidity axis is populated with fireballs following 
a gaussian distribution function given by p(i/fb) where i/fb is the rapidity of the fireball. 
Particles will appear when the fireball freezes out and will follow a thermal distribution 
centered around the position of the fireball 



The momentum distribution of hadron i is then calculated from the distribution of 
fireballs as given by Eq. [1] along the rapidity axis as follows 



where Ei—p^ is the the distribution of hadrons from a single fireball. The temperature 
T and the baryon chemical potential fiB will depend on the rapidity of the fireball and 
are not assumed to be constant. 

An important parameter is the width of the distribution. For the RHIC data at 
200 GeV this was determined from the 7r + 's as these are very sensitive to the value of 
a and less to variations in the baryon chemical potential. The width of the distribution 
a = 2.25 is compatible with the values quoted by the BRAHMS collaboration! [6J, 
e.g. a n + = 2.25 ± 0.02 and a n - = 2.29 ± 0.02. The hadrons described by eq. [2] are 
mainly resonances. Only a fraction of these are stable under strong interactions. The 
majority of them will decay into stable hadrons at chemical freeze-out, hence the need 
to implement multi-particle decays. 




(1) 




(2) 



d 3 Ni 
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3. Freeze-Out Curve 

We assume that the temperature T and the chemical potential \xb are always related via 
the freeze-out curve as given in Fig. [T] ; if the temperature varies along the rapidity axis, 
then also the chemical potential will vary. Thus a decrease in the temperature of the 
fireball will be accompanied by an increase in the baryon chemical potential. In other 
words, we assue a universality of the chemical freeze-out condition. This relationship 
between temperature and baryon chemical potential is very reasonable since all particle 
abundancies measured so far follow it. Once the width of the distribution of fireballs 
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Figure 2. Fit to the pion distribution as measured by the BRAHMS collaboration. 

has been fixed, we can go on with the dependence of the baryon chemical potential on 
the rapidity of the fireball (units are in GeV). 

fj. B = 0.026 + a y 2 FB (3) 

The dependence of the p/p ratio on the parameter a is shown in Figure [3]. The 
best value is a = 0.012 GeV. 

The variation of the temperature along the rapidity axis is shown in Fig. HI 
The temperature is maximal at mid-rapidity and gradually decreases towards higher 
(absolute) values of the rapidity. The baryon chemical potential follows a different 
pattern: it has a minimum at mid-rapidity and increases quite substantially towards 
higher values of the rapidity. Of particular interest are the largest values of the rapidity 
(say y ~ 4). It can be seen that there is a region of overlap with the data obtained at 
the NA49 experiment at CERN-SPS. 

4. Summary 

Particle yields measured in heavy ion collisions show an overwhelming evidence for 
chemical equilibrium at all beam energies. The rapidity dependence of the thermal 
parameters T and can now be determined over a wide range of rapidities and show 
a systematic behavior towards an increase in \xb as the rapidity is increased. 
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Figure 3. The p/p ratio as a function of rapidity. The curves show the dependence 
on the baryon chemical potential on the rapidity as discussed in the text. 




Figure 4. Values of the chemical freeze-out temperature and baryon chemical 
potential as a function of rapidity at the highest SPS (dotted lines) and RHIC (full 
lines) energies. 
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